BISPECTRAL COMMUTING DIFFERENCE OPERATORS FOR 
MULTIVARIABLE ASKEY- WILSON POLYNOMIALS 



PLAMEN ILIEV 



Abstract. We construct a commutative algebra Az, generated by d alge- 
braically independent g-difference operators acting on variables 21, Z2, . . . , z^i 
which is diagonalized by the multivariable Askey- Wilson polynomials Pn{z) 
considered by Gasper and Rahman [6]. Iterating Sears' 403 transformation 
formula, we show that the polynomials Pn{z) possess a certain duality be- 
tween z and n. Analytic continuation allows us to obtain another commuta- 
tive algebra An , generated by d algebraically independent difference operators 
acting on the discrete variables ni , n2, . . . , n^j, which is also diagonalized by 
P„{z). This leads to a multivariable g-Askey-scheme of bispectral orthogonal 
polynomials which parallels the theory of symmetric functions. 



Contents 

1. Introduction 1 

2. Multivariable Askey- Wilson polynomials 3 

2.1. Notations and one-dimensional theory 3 

2.2. Multivariable extensions 4 

3. g-DifFerence operators in 6 

3.1. Basic definitions 6 

3.2. The triangular operator Cd 7 

4. The commutative algebra Az 9 

4.1. Self-adjointness of Cd IC 

4.2. Construction of Az 15 

5. Bispectrality 17 

5.1. Duality of the multivariable Askey- Wilson polynomials 17 

5.2. The commutative algebra An IS 
Acknowledgments 21 
References 21 



1. Introduction 

In [1] Askey and Wilson introduced a remarkable family {pn{x) : n E Nq} of 
orthogonal polynomials on the real line depending on four parameters that satisfy 
a second-order g-difference equation 

A{z)Eg^z{Pn{x)) + B{z)pn{x) + C{z)E-l{pn{x)) = A„p„(a;), (1.1) 
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where x = ^(^ + 7)1 Eq,z is the g-shift operator acing on functions of z by 
= g{zq^^), v4(z), B{z), C{z) are independent of n, and A„ is inde- 
pendent of z. The orthogonahty iniphes that the polynomials Pn{x) satisfy also a 
three-term recursion relation 

anPn+l{x) + bnPn{x) + C„p„„i(x) = XPnix), (1.2) 

where a„, 6„ and c„ are independent of x and cq = 0. 

If we set = A{z)Eq^z + B{z)ld + C{z)E~\ then the operator acts naturally 
on the vector space consisting of all polynomials in the variable x and equation 
means that the polynomials Pnix) are eigenvectors of this operator with eigenvalues 
An. Similarly, let En denote the shift operator acting on an arbitrary function /„ = 
f{n) of a discrete variable n by En(fn) = fn+i and let L" = a„i5„ -I- 6„Id -I- CnE~^. 
Then the operator L" acts on the vector space of complex-valued functions of a 
discrete variable n S No and equation (|1.2[) means that Pn{x) are eigenvectors 
of the operator with eigenvalue x. Following [3] we can say that the Askey- 
Wilson polynomials solve a q- difference- difference bispectral problem. Moreover, 
this bispectral property can be used to characterize the Askey- Wilson polynomials 

pnni- 

As the title of the paper [1] states, the Askey- Wilson polynomials generalize 
those of Jacobi. In fact most of the polynomials used in analysis and mathematical 
physics can be obtained as special or limiting cases of these polynomials, see for 
instance [14] . For that reason, the Askey- Wilson polynomials are at the very top of 
the g- Askey scheme. For the terminating branch of the Askey scheme, the bispectral 
property has a natural interpretation within the framework of Leonard pairs [20j . 

An interesting question arising from the above discussion is whether there is a 
multivariable g- Askey scheme and if the polynomials in this scheme possess the 
bispectral property. We refer the reader to the book [4] for a very nice account 
of the general theory of orthogonal polynomials of several variables, which also 
indicates the two possible ways to proceed. 

One possible extension is linked to the theory of symmetric functions. In this set- 
ting, polynomials which are eigenfunctions of g-difference operators were proposed 
by Macdonald p/7j and Koornwinder [15] , see also [22] for connections with the bis- 
pectral problem mentioned above and [23] for multivariable g-Racah polynomials. 
These polynomial systems are associated with root systems and the corresponding 
polynomials are invariant under the action of the Weyl group. 

In a different vein, we can look for multivariable extensions, within the usual 
theory of polynomials of several variables, i.e. we consider an inner product in 
R'' and we apply the Gram-Schmidt process to all monomials, respecting the to- 
tal degree ordering. Within this context, multivariable Askey- Wilson polynomials 
{Pd{n; x;a) : n e Nq, x G M''} depending on d -I- 3 parameters ao, a2, ■ ■ ■ , a.d+2 
were discovered by Gasper and Rahman [6] . These polynomials represent g-analogs 
of the multivariable Wilson polynomials defined by Tratnik [21] . 

In the present paper we prove that the polynomials of Gasper and Rahman are 
also bispectral. More precisely, we define two commutative algebras Az and An 
of operators which are simultaneously diagonalized by the polynomials Pd{n] x; a). 
The algebra Az is generated by d algebraically independent g-difference operators 
£f , . . . , £^ acting on the variables zi, Z2, ■ ■ ■ , Zd where Xj = ^{zj + z~^). Thus the 
operators £| act naturally on the vector space C[xi, . . . ,Xd] and are diagonalized by 
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the polynomials Pdin^x^a). Likewise, the algebra An is generated by d algebraically 
independent difference operators £",..., £2 acting on the variables ni, . . . ,nd and 
the polynomials Pd{n;x;a) are eigenfunctions of these operators (considered as 
elements of the vector space of complex- valued functions defined on Nq). 

The paper extends the results in our joint work with Geronimo 8J and it raises a 
series of interesting questions. For instance in the one-dimensional case, for specific 
values of the free parameters, the Askey- Wilson recurrence operator L" commutes 
with a difference operator of odd order. Jointly with Haine [H] we explored this 
fact to connect the above theory to algebraic curves, having specific singularities, 
or equivalently, to specific soliton solutions of the Toda lattice hierarchy, using 
the correspondence established by van Moerbeke-Mumford [HI US] and Krichever 
[16j . In particular, this led to a different explanation of the bispectral property 
using algebro-geometric considerations [11]. Moreover, this approach suggested us- 
ing techniques from integrable systems (such as the Darboux transformation) to 
construct extensions of the Askey- Wilson polynomials which satisfy higher-order 
g-difference equations |12j . In the multivariable case, algebro-geometric methods 
were used by Chalykh [2], within the context of symmetric functions, to give more 
elementary proofs of several of Macdonald's conjectures. It is a challenging prob- 
lem to construct a Baker-Akhiezer type function for the operators discussed in 
this paper and prove the duality using algebro-geometric tools. Another interest- 
ing question is to classify the multivariable orthogonal polynomials that satisfy 
g-difference equations. Recently discrete multivariable orthogonal polynomials sat- 
isfying second-order difference equations were classified in our joint work with Xu 

m- 

The paper is organized as follows. In Section|2]we recall the basic definitions and 
orthogonal properties of the one-dimensional Askey- Wilson polynomials as well as 
the multivariable extension proposed by Gasper and Rahman. We also show that 
a change of variables leads to the multivariable g-Racah polynomials discussed 
in [7] . In Section [3] we define a g-difference operator Cd acting on the variables 
zi, Z2, ■ ■ ■ , Zd- We prove that it preserves the ring of polynomials in xi, X2, . . . ,Xd 
and we establish its triangular structureQ. In Section |4] we show that Cd is self- 
adjoint with respect to the inner product defined by the multivariable Askey- Wilson 
measure. This leads to the construction of the commutative algebra Az diagonalized 
by the multivariable Askey- Wilson polynomials. In Section [5] we iterate an iden- 
tity of Sears to show that, appropriately normalized, the polynomials Pd{n]x;a) 
possess a certain duality between the discrete variables (ni, n2, . . . , n^) and the 
continuous variables (zi, Z2, . . . , Zd)- Using this duality we obtain the commutative 
algebra An of difference operators acting on the variables ni, 712, . . . , which is 
also diagonalized by Pd{n]x;a), thus proving the bispectrality. 

2. Multivariable Askey- Wilson polynomials 

2.1. Notations and one-dimensional theory. Throughout the paper we assume 
that (7 is a real number in the open interval (0, 1). We shall use the standard q- 
shifted factorials 

n — 1 00 k 

{a;q)n = ]^(l-ag''), {a;q)oo = JJ(l-ag'') and (01,02, . . ■,ak;q)n = (%;?)«, 

k=0 k=0 j=l 



with respect to the total degree ordering 
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and the corresponding basic hypergeometric series, see j5j for more details. We 
denote by Pn{x] a, 6, c, d) the Askey- Wilson polynomials [1] 



{ab,ac,ad;q)n , 
p„(x; a, 6, c, d) = 403 



q ^.abcdq^ ^,az,az ^ 
ab, ac, ad 



(2.1) 



where x = \{z + i). When a,b,c,d are such that max(|a|, |c|, < 1 the 
polynomials defined by (|2.ip satisfy the orthogonality relation 

1 /"^ , , ,N / ,\W(z;a,b,c,d) c i ^\ 

p„{x;a,b,c,d)pjn[x;a,b,c,d) ^ — dx = dn^mK, (2.2) 



where 



and 



27r J_i VI - a: 



(az,az '^,bz,bz '-,cz,cz ^,dz,dz ,q)oo 
^ ^ {abcdq''-^;q)n{abcdq^''; g)oo 4) 



2.2. Multivariable extensions. Consider x = {xi,X2, ■ ■ ■ ,Xd) G M.'^. For jsjl < 1 
we put 

xj=cos(0,) = ^(z, +Z-1), (2.5) 

where Zj = e^^' for j ~ 1,2, ... ,d. We define a weight function w(z) = Wd{z; a) 
depending on c? + 3 nonzero real parameters ao, ai, . . . , ad+2 by 

-(^) = "^-^^^^'^^'^^^7 , (2.6) 

with the convention that zq — ao and 2^+1 = ad+2- We shall assume that the 
parameters ak are such that if |zj | = 1 for j = 1,2, ... ,d then 

Ittfe+iajT^z^^^z^^l < 1 for k = 0,l,...,d. 

It is easy to see that this is equivalent to the constraints 

< \ad+i\ < \ad\ < • • • < < min(l, \ao\^) 



\ad+i\ ^1 , \ad\ 

-. r- < \ad+2\ < -. r 

\ad\ \ctd+i\ 

Let us consider the following measure on [—1, 1]° 



(2.7) 



and the corresponding inner product 



(2^) f=\ Jl - a;? 



if, 9) = / f{x)gix)d^l{x). (2.8b) 
J[-i,i]'' 

The following theorem established in [6] constructs an explicit basis of polyno- 
mials mutually orthogonal with respect to the inner product defined by p.6p - (|2.8p . 
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Theorem 2.1. For n = (ni, n2, ■ ■ ■ , rid) G Nq define 

Pd{ri;x;a) = f[p^^ f:,^■a,q'''~^,^q^'~^,^z,+u^z~l^ , (2.9) 
jJi \ "o / 

where Nk ^ ni + n2 + ■ ■ ■ + and Nq = 0. // the parameters {aj}'^^Q satisfy (12. 7p . 
t/ien for every n, m £ Nq we /laue 

{Pd{n;x;a),Pd{m;x;a)) ^ Sn.niHn, (2-10) 

where the inner product is defined by ()2.6p anrf p.8|l . and 



(2.11) 



-t^n — 11/ 2 2 



X 



E{ 1,1} ("'i+ia|+2'Z^",a<i+iad+2«o 9 '';9)oo 
For the convenience of the reader we sketch the proof. 
Proof. We can write the weight (|2.6p as w(2;) = wi (2)71)2(2:), where 

wi(z) = — J — , (2.12a) 

is the one-dimensional Askey- Wilson weight for zi, depending also on 22 and 

W2{z) = — n,t2(^|.^7^;g)oo ^2.i2b) 

nfc=2nei,e2G{-l,l} ("^+1"^ 4+A''^l)^ 

is independent of zi. Notice also that in the right-hand side of equation (|2.9p only 
the first polynomial Pm depends on 21. Writing the integral representation of the 
inner-product in (j2.10p we see that the only terms depending on xi are Pm , Pmi 
and wi{z)/ y'^^^^xf. Using (|2.2p we obtain 

1 w,{z) ^ u ~ r ^ 

TT / PniP7ni—? „ dXi = 0ni,miftniWi(22), 



where 
and 



ri, \"o 



Wl(z2) = Jl 



se{-i,i} (^a2g"i2|,^g»i2f;g 

It easy to see that wi(22)u'2(2) coincides with the weight Wd-i{z,a) given in (|2.6p 
for the variables 21 = 22,- • • j-2d-i — Zd and parameters fio = ao, fii = Q;2<z"^, 
<52 = a39"S- • • 7 <5d = oid+iq"'^ , otd+i — ad+2- Moreover, if we denote n = 
(n2, na, . . . , Ud) € Nq^^ one can check that Pd-i{n\ i; a) coincides with the product 
0^=2 P"3 consisting of the last {d— 1) polynomials for Pd{n\x\a). The proof now 
follows by induction. □ 
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Remark 2.2. One can easily obtain the parametrization used by Gasper and 
Rahman [6], by introducing parameters a,b,c,d, a2,a3,...,ad related to aj by 
the formulas a = ai, & = cti/al, c = ad+iad+2/ad, d = aa+i/ {ada.d+2) and 
a/c+i = ak+i/ak for fc = 1, 2, . . . , d - 1. 

Remark 2.3. The multivariable g-Racah polynomials discussed in [7j can be ob- 
tained as follows. Suppose 

^ = g^, where N 
ad+i 

Making the substitution 

Zk = akq^'' for fc = 0, 1, . . . , d + 1 
we see that polynomials in ()2.9I) become the multivariable g-Racah polynomials 

Rn = \{ (vk - iV,_i; 4g'^^-^-\ a^^g^-^+^-S y.+i - 

(2.13) 

where are the one dimensional g-Racah polynomials defined by 



rfc(?/; a, 6, c, A^) ^ {aq, bcq, q ^ ]q)k{q^ / cf'"^ a<Pz 



'q-'',abq''+'^,q-y,cqy-^ 

L -N '^Q^Q 
aq,bcq,q 



The polynomials i?„ are orthogonal on = yo < yi < 1/2 yd l£ Vd+i — N 

with respect to the weight 

/ N YJ i(^l+l/o^hl)vk + l-Vkio^k+l'-'l)vk+l+yk -rr,, 2 2yk\ f '^l<:~l\^^'' 

P{y)=\\ 7 ^ 7 — 2-^^ II l-afc?*^) 

[l-^QJvk+i-ykiQo^k^Qhk+i+Vk j^Jl \ o^k J 

The parametrization in [7] can be obtained if we replace ao, ■ ■ ■ , ctd+i by ai, . . . , 0^+1, b 
related via the formulas ai — af, Uk — ct'k/'^'k-i ^'^^ ^ ~ 2, . . . ,d + 1 and b = 
a? /(gag). 

Thus, all difference equations derived for the multivariable Askey- Wilson polyno- 
mials later in the paper, can be translated to the multivariable g-Racah polynomials, 
using the change of variables and parameters above. 

3. g-DlFFERENCE OPERATORS IN 

In this section we construct a g-difference operator Cd which is triangular and 
self-adjoint with respect to the inner product 



3.1. Basic definitions. Consider the ring Vz±^ = C[zf^, zf^, . . . , z^^] of Laurent 
polynomials in the variables zi, Z2, . . . , Zd with complex coefficients. We denote by 
Vx — C[a;i, X2, ■ ■ ■ , Xd] the subring of 'Pjii consisting of polynomials in the variables 
xi,X2, ■ ■ ■ ,Xd, where 

X, = i(z, + Z7I) for i = 1, 2, . . . , d. (3.1) 
For J = 1, 2, . . . , d we define an automorphism Ij of Vz±i by 

Ij{zj) = z^^ and Ij{zk) = Zk for j ^ k. (3.2) 

Clearly, Ij is an involution (i.e. Ij o Ij = Id) which preserves Vx- Conversely, if 
a polynomial p G Pz*! is preserved by the involutions Ij for j — 1, 2, . . . , d then 
V&Vx. 



BISPECTRAL OPERATORS FOR MULTIVARIABLE ASKEY-WILSON POLYNOMIALS 7 



Let {ei, 62, . . . , Crf} be the standard basis for C^. We denote by Eq^ z. , Aq, z. and 
Vq,z^, respectively, the g-shift, forward and backward difference operators in the 
j-th coordinate acting on functions f{z) as follows 

Eg, zj{z) = /(zi, Z2, . . . , Zj-g, . . . , Zd) 

Aq.zJ{z) = {Eq.z^-l)f{z) 

Vq.,zjiz) = {l-E-X)fiz). 
Throughout the paper we use the standard multi-index notation. For instance, if 



(i^i, . . . , fp) G Z'' then 



d' ^q.z—-^q.Zi-'^q,Z2 -'^q,z,,-, ZQ - (ZlQ , Z2q , . . . , ZdQ ) 



1^2 ' ' ' •^d ' ^q, z ^q, zi ^q, Z2 ' ' ' ^q, z^ 

and \v\ = vi + V2 ^ h i^d- 

We denote by ~ C(zi, . . . , Zd)[£'^^^ , . . . , iS^^^] associative algebra of 
g-difference operators with rational in z coefficients, i.e. the elements of I?z are 
operators L of the form 

L = Y,Uz)El,, 

where S* is a finite subset of Z'' and lu{z) are rational functions of z. Thus, the alge- 
bra Vz is generated by rational functions of z, the shift operators Eq, zi , Eq, z^ , ■ ■ • , Eq, zd 
and their inverses E~^\^ , E~\^ , . . . , E~^^^ subject to the relations 

Eq,z, ■g{z)^g{zq'^)Eq,z^, (3.3) 

for all rational functions g{z) and for j ~ \,2, . . . ,d. For every fc S {1, 2, . . . , rf} the 
involution Ik can be naturally extended to I?z, by defining 

h {Eq., z, ) = E-^ and h {Eq. z, ) = Eq^ z, for 3 + ^- (3.4) 

Indeed, it is easy to see that Ik is correctly defined because the relations p.3p are 
preserved under the action of /fc, i.e. we have 

lk{Eq,z,) ■ h{9{z)) = 4(5(zq^0)4(i?,,z,), 
for e {l,2,...,d}. 

For the sake of brevity we say that i £ I?z is /-invariant if Ij [L] = L for all 
jG {l,2,...,d}. 

For fc e No we denote by Vj^ the space of polynomials in Vx of (total) degree at 
most k in the variables xi, X2, ■ ■ ■ , Xd, with the convention that = {0}. 

Definition 3.1. We say that a linear operator L on Vx is triangular if for every 
e No there is Cfc G C such that 

L{p)=CkP mod T'^-i for all p e 7^^ 

3.2. The triangular operator Cd- Below we define an operator which is I- 
invariant and triangular. Clearly, an operator which is /-invariant is uniquely deter- 
mined by its coefficient of E^ ^ (or equivalently A'^ ^) with I'i > for j = 1, 2, . . . , d. 

Let — (i^i, z^2, ■ • ■ , i^d) € {0, {0}'', and let {i^i^, Vi^, . . . , i^i^} be the nonzero 
components of v with 1 < ii < 12 < ■ • ■ < is < and s > 1. Denote 

Ai, =(1 - ai^Zi^){l - ai^ZiJal) 

X ^1 ITTTi ITTl l^o.oaj 



m=i(i-4)(i-94) 

X (1 - ad+iad+2ZiJaiJ{l - ad+iZiJ{ai^ad+2))- 
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An arbitrary u ^ can be decomposed as v = — ^ where e Ng with 

components z/^ — max(i'j, 0) and vj = — min(j/j,0). For v e { — 1,0, {0,1}'' 
we define 

A, = r"(A,++,-). (3.5b) 

Here I'^ is the composition of the involutions corresponding to the positive coor- 
dinates oi v~ . 

FinaUy, we define the operator 

Cd^ Cd{zi,...,Zd;ao,ai,...,ad+2)= ^ (-1)''^ '^i^'^q.^z'^g, z- 

i^e{-i,o,i}''\{o}'' 

(3.6) 

Since Ii{Aq^zi) = — Vq,z. and Ii{Aq,zj) ~ ^q.,zj for i ^ j, the operator defined 
above is /-invariant. 

Lemma 3.2. Let L G be an I -invariant q-dijference operator. If 
d 

\{{l-z^^)L{p)eV,±^ for every p^Vx^ (3.7) 

then L preserves Vx, i-e. L{Vx) Vx- In particular, the operator Cd defined by 
p.5p - p.6p preserves Vx- 

Proof. Let p{x) £ Vx- Since both L and p are /-invariant, it follows that L{p) is 
also /-invariant and therefore to show that L{p) S Vx it is enough to prove that 
Liji) E Vz±i- From (|3.7p it follows that we can write L{p) as a ratio 

L{p) = ^ , where pi(z) G V^±i. 

Notice that Ij{zj — z~^) = —{zj — zj^), and therefore Ij{pi{z)) = —pi{z) for every 
j = 1, 2, . . . , d. This implies that pi{z) = when Zj = ±1 and therefore pi{z) is 

divisible (in the ring Vz±i) by nj=i(^i ~ -^7^)' Pi'^'^i^S that L{p) G Vz±i. 

It remains to show that Cd satisfies the conditions of the Lemma. This follows 
easily from formulas (|3.5p - (|3.6p combined with the fact that Aq^zjix'j) is divisible 

by Aq^zAx,) = %^7?^: and Vq^z^ix'^) is divisible by V,,z,(x,) = 

for every fc G N. □ 



Proposition 3.3. The operator Cd defined by p.5p - p.6p is triangular. More pre- 
cisely, we have 

Cd{p{x)) = -(1 - q-^) (^1 - ^^9''"') Pix) mod Vt^ for every p{x) G V^ 

(3.8) 

Proof. From Lemma [3.21 we know that Cd{p{x)) G Vx for every p{x) G T^^;. It is 
enough to prove that (|3.8p holds when p(a;) — x"" — a;"^ ^ • • • x^'' . Clearly the 
highest (total) power of x in Cd{x") can be uniquely determined by the highest 
power of z. Thus it is enough to see that, after canceling the denominators of the 
coefficients A^, of the operator Cd, we have 

Cd{xn = -^(1 - g-l«l) (l ^gl-l-^ 



a linear combination of z'' with < 
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Let 7^ e {-1, 0, 1}''. Notice that 

= ^(9"^- - 1)^? + Oiz;^-') (3.9a) 
= ^(1 - q-"^)^' + 0{z;^-') (3.9b) 

and therefore A[^+ Vg^zl^") is a hnear combinatfon of z" and z''' where \k\ < \n\. 
From (|3.5p it follows that the denominator of A^, is a polynomial of degree 4s in 
zi, Z2, ■ . . , Zd^ where s denotes the number of nonzero components of v. However, 
the numerator is of degree 4s if and only if all the nonzero coordinates of v have 
the same sign (i.e. they are all positive, or all negative). If at least two of the 
coordinates of i' have different signs, then the degree of the numerator is less than 4s 
in the variables zi, Z2, . . . , Zd- Thus we conclude that Cdix"^) is a linear combination 
of 2" and z'' where |fc| < |n|. Moreover, the coefficient of z" can be computed 
by extracting the coefficient of z" from the terms (—1)'"^ 'A^A^^^VJ^ zi^") 
ly G {0,1}'^ U {—1,0}'^ in the representation (|3.6p of the operator Cd- Using 
formulas (|3.9p we deduce that 

Cd{x^) = ^j^-^" + hnear combination of z'^ with < 

where Cn — + with 

4 = ^ E n (^"'-D (3-lOa) 

such that i^s — 1 

c-= E n !)■ (3-lOb) 

0=£i^e{0,~l}'' sG{l,...,d} 

such that i^s — — 1 

It is easy to see (for instance by induction on d) that 
05^i/g{o,i}'i se{i,...,<i} 

such that i^s — 1 

which combined with p.lOp gives 

Q,2 

c+ = ^(gl"l-l) and c-=q-H-l. 



Thus 



c„---(l-<Z-l"l)(l-^^l 



completing the proof. □ 



4. The commutative algebra Az 

In this section we show that the operator Cd defined in Section [3] is self-adjoint 
with respect to the inner product (|2.6p - (|2.8p . This allows us to construct a commu- 
tative subalgebra Az of Vz , generated by d algebraically independent g-difference 
operators, which is diagonalized by the polynomials (|2.9p . 
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4.1. Self-adjointness of Ld- We start with a simple lemma. 

Lemma 4.1. Consider an inner product defined by equations ()2.8|) and let L — 
"Yliu^s 2 G U for every v ^ S the following conditions are satisfied 

(i) Cu{z)w{z) = C^,{zq'')w{zq''); 

(ii) The function C^{z)w{z) is holomorphic on the d- dimensional torus 

= {z G C : \zj \ = q-^''^ for j = 1, 2, . . . , d} for every t G [0, 1] 
then L is self-adjoint with respect to the inner product ()2.8|) . 

Proof. Using the change of variables (12.511 we can write the inner product as an 
integral over the d-dimensional torus T = {z e : |zj | = 1 for j = 1, 2, . . . , d} 

^/'^'^ = 71^ / ri — • (4-1) 

{Amy J J ^J-^ Zj 

Denote = C,y{z)E^^^. Then from (|4T|l we obtain 

(47ri)''7T -^^1 Zj 

replacing z by zg" 

r dr- 



{Am) 



using (i) 



1 



using (ii) we can replace T!j' by Tq = T 

= 71^ / C.(z)/(zg^)g(zWz)n^ 
= {L.f.g). 

The proof follows immediately by writing L as a sum of L^^s. □ 



In order to apply Lemma 14.11 we need to rewrite the operator Ca as a linear 
combination of the q-shift operators i?^ ^ . 

For j e {0, 1, . . . , d} and (fc, e {0, l}^ we define as follows 

„2 



i?f (z) = 1 + - ^7^:7;^"^ (4.2a) 

^0,1^^^ ^ / _ «.-+i^.-^.+A ^ _ ^^±i£2±i\ (4.2b) 



j+i 4aj+ia;ja;j+i 



a," 



1,0^^^ ^ / _ a,+iz^\ / ^ _ 2c) 
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In the above formulas Xj is related to Zj by p.Sp . with the convention zq = ao and 

Zd+i = Cid+2- We extend the definition of B^'^{z) for (fc, I) e {—1, 0, 1}^ by defining 

BJ^\z) = Ij{bY{z)) for ; = 0, 1 (4.2e) 

Bf-\z) = I,+,{Bf\z)) for fc = 0, 1 (4.2f) 

Bj'^-\z)^I,{I,+,{B]^\z))). (4.2g) 

Next, for j e {1, . . . , d} we denote 

6°(z) = (l-(7^|)(l--?z-2) (4.3a) 

b]{z) = {l^zj){l-qz^^) (4.3b) 

6-i(z)-/,(&](z)). (4.3c) 

Finally, for u e {-1, 0, l}"* we put 

a(z) = (,(g + l))'^-|^"l-l-l ; \ (4.4) 
with the convention i^q = t'd+i = 0. 

Proposition 4.2. T/ie operator Cd — Cd{z;a) defined by p.5p - (|3.6p can 6e written 
as 



(4.5) 



where C^{z) are given by (|4.2p . (|4.3p and (|4.4p . 

Proof. We shall prove the statement by induction on d. For d — 1 formulas p.5 
give 

(l-z2)(l-gz2) 

(1-Z2)(5_^2) 



= /l(^ei) - 

Using (112]), (1131) and ((13) we get 

(l-zf)(l-gz?) 

C_ei =/l(Cei(z)) - A-ei 
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where xq — ^(cto + cto = + ^i ^) and X2 — \{az + ^). We need to 

check that 

C,=A,,{E,,,,-1)-A_,,{l-E-X) 

V ic^o {q + l)aoJ 

which amounts to checking that 

4a22:o2:2 



-Ae, - A_ei = Co - 1 + ^ 

This equahty can be verified by a straightforward computation using the exphcit 
formulas for A^-^ , A-^i and Co above. 

Let now d > 1 and assume that the statement is true for d—1. We can write Cd 
as foUows 

( _ ad+iad+2 \ f ^ "d+l A 

I f-ll \ "dZd J \ adad+2ZdJ I ^ s „,„ 

where £", are g-diff'erence operators in the variables zi, Z2, . . . , Zd-i with 
coefficients depending on zi, . . . , Zd and the parameters aoi cii, . . . , q;c(+2- Clearly, 
the operators are uniquely determined from Cd- This imphes that 

they are /i,/2, . . . ,Id-i invariant and therefore they are characterized by the co- 
efiicients of ^ = ^qUi^q%2 ' ' ' -^qX-n where z = (zi, Z2, • ■ • , z^-i) and P = 
(j^i, 1^2, ■ ■ ■ , '^d-i) G {0, 1}''^^. The coefficient of A° j in £' is equal to {l~adZd){l~ 
adZd/ctQ) (it comes from the term A^^Aq^zd in ^d)- Using equations p.Sp one can 
see that the coefficients of j for j> G {0, 1}''~"'^\{0} are the same as the coefficients 
of the operator Cd~i{zi, . . . , Zd-i; c^o, ■ • ■ j cid-i, ct'^, o^'d+i) where a'^ — adZdq^^^ and 
^'d+i = q^^^^- Notice that 



qal (q'+l)Q;o 



(1 - adZd)(l - adZd/a^) 



Thus, using the induction hypothesis, we deduce that 

C ^Cd-i{zi, . . . , Zd-i; ao, . . . , ad-i, adZdq^^'^,q~^^'^) 
+ (1 - adZd)il - adZd/al) 

~ CvEq,z^ 
i/G{-l,0,l}''-i 



(4.7a) 



where C^ are computed from (|4.2|) - (|4.4[) for the operator Cd-i with parameters 
given in (|4.7ap . For C" we have 

C"^Id{C')^ Id{Cl)El,^ (4.7b) 

!/£{-!, 0, 1}''-! i/e{-i,o,i}''-i 
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For the last operator we obtain 

= £d~i{zi, . . . , Zd-i;ao, . . . , a^-i, a^+i, ad+2) 

V C'E"" - (1+ "^+1 - "^^d+iXqXd+A (4.7c) 
i^e{-i,o,i}'*-i ° ^ " / 

where C"' are computed from ()4.2p - (l4.4p for the operator £d-i with parameters 
given in (|4.7c|) . 

Notice that the last term in ()4.7c|) gives the last term in the right-hand side of (14. 5|) . 
Thus, it remains to show that we get the stated formulas (|4.4p for the coefficients 
Cj, in (|4.5p . using the decomposition (|4.6p . Since the operator Cd is /-invariant, it 
is enough to prove the formulas for Cjy when v has nonnegative coordinates. We 
have two possible cases depending on whether Vd = ov Vd = 1. 

Case 1: = 1- Write v = {v' , 1) with v' e {0, 1}'^"^ From (j46|l it is clear that 
E'^ ^ appears only in the first term on the right-hand side and we have 



ad 



{l-zl){l-qzl) 



Notice that the factors (^B^^*"'^''+^y jn formula (|4.4p for C^, are the same as the 
factors B^f^'""''^^ in formula (|4.4p for d, when k = 0, 1, . . . , d - 2. Similarly, the 
factors {b^i!')' in formula (|4.4p for C^, are the same as the factors fe^*" in formula 
(|4.4p for Cu when k = 1. This combined with the last formula above and 

h\ = {i- ,4). BY =(i- '^^±1^^) ( 1 _ i!^) 

V Ud J \ OidOid+2) 

show that in order to complete the proof we need to check that (B"^^^'^y — B"^^^''^ . 
This can be easily verified from the defining relations (14. 2p by considering the two 
possible cases Vd-\ = and i^d-i = 1- 

Case 2: Vd = 0. Let us write again v = (i^',0) with v' e {0, 1}''^^. Then from 
(|4.6p we deduce 

/ 2 „ £d+l£d+2£dA ( 1 _ ad+l^d \ 



(l-zj)(l-(7^^) 



rill V "dZd y V adad+2Zdy „,„ 

We need to check formula (|4.4p for C^. Again the factors Sj^^'' '"''^^ for fc = 0, 1, . . . , d- 
2 and the denominator factors l)'^ for fc = 1, 2, . . . , c? — 1 are common for Cv, C^/, 
C'li and C",'. Thus we need to verify that 



B 



d-\ _ /n'^d-l.Oy V / V aaaa + 2j 



''^''^'\l-qz])il-qz-^)- (1--?)(1-'Z^^) 



1 - 



Qd+lQd+2 



fl 2£l±^') 

^ ni^d-l.Ow/ V "t'^'' J\ adad+2Zd J / Tj'^d-lfi\lll 

y^d-i ) — ^277; ^ y^d-i ) ■ 
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Using the explicit formulas ()4.2|) and considering separately the two possible cases 
i^d-i = and lyd-i = 1 one can check that the above equality holds, thus completing 
the proof. □ 



Proposition 4.3. The operator Cd defined by ()3.5|) - (I3.6|) is self-adjoint with respect 
to the inner product (|2.6p - (|2.8p . 



Proof. We shall use Proposition 14. 21 to check that the conditions in Lemma HTT] are 
satisfied. We can represent the weight w{z) in (|2.6p as 

where 

w'kiz) ^ {zl,z~^;q)oo (4.8a) 

Wk(z)^{ Zk+iZk, , , ;g . 4.8b) 

V Oik Oik Zk Oik Zk+i ak ZkZk+i ) ^ 

Using Proposition 14. 21 and formula (14. 4p for the coefiicients of the operator Ld{z\ ol) 
we see that condition (i) of Lemma 14.11 will be satisfied if we can show that for 
every v E { — 1, 0, 1}'' we have 



(4.9a) 



hl'^'izq'^) w'^{zq-) 
and 

- «=> ,4.9b) 



It is easy to see that if these formulas are true for some v then they are also true 
for —V. Since b^f!'{z) depends only on Zk, it is enough to check (|4.9ap for Vk = 
and Vk = 1- The case = is trivial because clearly both sides of equation (|4.9a[) 
are equal to 1. If = 1, then the right-hand side of (|4.9a|) is 

<{z) ^ {zl,z-^-q)^ ^ {l^zl){l~qzl) 
Kizq'^) (q'zl z-\-^; q)^ f 1 _ _ 

Using (j4.3bp - (|4.3c[) we see at once that the last expression equals b\{z) /b^^ {zq'^'') 
completing the proof of (|4.9ap . Similarly, B'^'"'^'''^^ (z) depends only on Zk and Zfc+i- 
Thus it is enough to verify (|4Jb| for (i^fc,i^fc+i) = {(1, 0), (0, 1), (1, 1), (1, -1)}, 
which can be done in the same manner by a straightforward verification using 
formulas ((i?^ . 

Condition (ii) in Lemma [Ol will follow if we show that for every v £ {—1, 0, l}"* 
the functions 

w'kiz) 



b-kHz) 

and 



for fc=l,2, (4.10a) 



for fc = 0,l,2,...,d (4.10b) 



<{z) 

are holomorphic on for < e [0, 1]. Notice that 



w 



',iz) = (1 - zl){l - qzl){l - z,'){l - qz-^){zlq\ z^^^- q) 
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which combined with (14. 3p shows that aU the factors in the denominator in (I4.10a|) 
cancel leading to a holomorphic function on T^. Next we show that function in 
(|4.10bp is holomorphic on TJ'. From (|4.9bp it follows that if this is true for j/, then 
it will also be true for —v. Thus we can consider only the cases: (vkT^k+i) = 
{(0,0), (1,0), (0,1), (1,1), (1,-1)}- The case (vk^Vk+i) = (0,0) is trivial since 
equations p.7p guarantee that w'l{z) has no zeros on the the d-dimensional torus 
T = {z S : |zj| = 1 for j — 1, 2, ... , d}. Let us consider for instance the case 
{vk, Vk+i) = (1, 1 ). Th en for z G we have 1 < \zkZk+i\ < and \zk/zk+i\ = 1. 
From (|2.7p and (I4.8b|) it follows that the only factors in w'l{z) that can vanish 
during the homotopy are 

^ _ ak+iZkZk+i \ A _ qak+iZk+iZk 



ak J \ cik 

From (j4.2dp we see that these two factors cancel with B]^^ showing that the function 
in (|4.10b[l is holomorphic on T^. The remaining three cases follow along the same 
lines. □ 

As an immediate corollary of Proposition 13. 31 and Proposition 14.31 we obtain the 
following theorem. 



Theorem 4.4. Let ao, ai, . . . , Q!c(+2 be real parameters satisfying (|2.7p and let 

{Q{n]X]a) : n G Nq} he a set of orthogonal polynomials with respect to the inner 
product (j2.8p . i.e. Q{n;x;a) is a polynomial of total degree \n\, orthogonal to all 
polynomials of degree at most |n| — 1. Then Q{n; x; a) is an eigenfunction of the op- 
erator Cd defined by p.5p - (|3.6p with eigenvalue ji = —(1 — (/^'"l) ^1 — -^T^f?'"'^^^ • 

4.2. Construction of Az- Next we focus on the polynomials Pd{n;x;a) defined 
by (I2T|) and (EH). 



Proposition 4.5. Let a G {C*y+^. For j G {1,2, . . . ,d} define 

£j = Cj{zi, . . . , Zj]ao, . . . ,aj+i, Zj+i) (4.11a) 



/ij = -(1 - + j I _ _J^qn,+n2 + -+n,-l j (4.11b) 



a, 







Then the polynomials Pd{n; x; a) defined by (12. ip and (|2.9p satisfy the spectral equa- 
tions 

£,jPd{n; x; a) = fJ-jPd{n; x; a) (4.12) 

for every j G {l,2,...,d} and the operators £j commute with each other, i.e. 
is a commutative subalgebra ofV^. 

Proof. If Qj are real satisfying (|2.7p , then (|4.12p for j = d follows immediately from 
Theorem 14.41 Notice that for fixed n G Nq both sides of (|4.12p depend rationally 
on the parameters aj . Thus when j = d equation (I4.12p holds for arbitrary values 
of the parameters aj . From (|2.9p it clear that for j < d the product of the first j 
terms on the right-hand side is precisely 

Pj{ni, . . . ,nj;xi, . . . ,Xj;ao, . • . , aj+i, z^+i). 



while the remaining terms Y\'i=j+iPni do not depend on the variables xi, . 
This shows that (|4.12p holds also for j < d. 



. X'j 
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It remains to prove that the operators £| commute with each other. Fix z 7^ j G 
{1,2,..., d}. From (|4l2l) it follows that 

[£f , £|]Pd(n; x; a) ^ 0, for all n e N^. 

This means that the operator L = e vanishes on Vx and we want to 

show that it is identically equal to zero, i.e. its coefficients are identically equal to 
zero. Equivalently, it is enough to show that for fixed k = (fci, . . . , kj,) £ Nq the 
operator L' = (which also vanishes on Vx) is identically equal to zero. We 

can assume that the components kj are large enough, and therefore the operator 
L' will contain only nonnegative powers of Eg, z- If we denote 

1 

then the operator L' can be uniquely written as 

L'= Uz)Dlz. 

i/|<M 

where M e N is large enough and lu{z) rational functions of z. Now we use the 
fact that L'{x"-) = for all n € Nq such that \n\ < M. We obtain a homogeneous 
linear system for l^{z) with determinant 

Det= det 

\i'\,\n\<M 

In the above determinant we order in the same way v (the rows) and n (the columns) 
respecting the total degree, i.e. v < jj, when |z^| < \^\. It remains to show that Det 
is not identically equal to 0. Notice that 

and therefore if we put 

_ n,ti(g"-'"-+^g).. 

then 

E>q.zX^ = Si^.nz""" + terms involving with |fc| < |n| - \u\. 
From this formula it follows that 

Det — det {gu.n) + a linear combination of z^ with |A:| < 0. 

\iy\,\n\<M 

Clearly, g^^n 7^ if and only if Vj < nj for all j. This shows that matrix {gu.n) is 
upper triangular with nonzero entries on the main diagonal, hence detu,n{gv,n) 7^ 
completing the proof. □ 
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5. BiSPECTRALITY 



5.1. Duality of the multivariable Askey- Wilson polynomials. Since g"^ de- 
termines rij uniquely modulo ) '^'^ consider in this section complex variables 



n = (rii, ...,nd) G (i 



mod 



j^^Zj , complex variables z = (zi, Z2, . . . , z^) e 
(C*)'^ with nonzero components, and nonzero parameters a — {ag , . . . , ad+2) G 

We define dual variables ri e mod j^^^ , ^ G (C*)'' and dual parameters 

€ {C*y+^ by 



Zi 



Ctd+l-jZd+l-j 
ad+2-jZd+2-j 

d+2 

an 



"^+2-J ^jVd+l-,-l/2 



1/2 

ao 



for J = 1,2, ... ,d 
for J = 1,2,..., (i 

for j = 1,2, . . . + 1 



where as before we have set Zd+i — Q-d+2 and = «i 
(|2.5|) we put = + z^^"^). 



(5.1a) 

(5.1b) 
(5.1c) 
(5.1d) 

(5.1e) 

rife . In analogy with 



Lemma 5.1. The mapping f : (n, z,Q!) — > (n, 5,5) given by (15. ip defines an invo- 



lution on ( C 



27ri 
log(9) 



X (C*) 



*\(i+3 



Proof. Using formulas (|5.ip one can easily check that f o f = Id. 

Iterating a formula of Sears we obtain the following lemma. 
Lemma 5.2. //fc e Nq and ahc — defq^^^ then 



□ 



{d,ej;q)k4<p3 



;q,q 



X 403 



= c'=(fe,agi ''/e,e/c;q)k 



(5.2) 



Proof. Sears' formula (see for instance 5, page 49, formula (2.10.4)]) gives 



{d,ej;q)ki(t)3 



;q,q 



X 403 



d,e,f 
q-^,a,d/b,d/c 
d, agi^Ve-a^'^V/ 



a'^Xd, e/aj/a,;q)k 



(5.3) 



Applying the same formula on the right-hand side of (j5.3p with a, 6, c, c?, e, / replaced 
by c?/&, a, d/c, og^^'^/e, d, aq^~^ / f respectively and using the identity 



(s;g)fc = (-s)'=9 



we obtain (|5.2p . 



□ 
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Let us normalize the multivariable Askey- Wilson polynomials p.9|) as follows 

P,(„; ^. a) = (ad+iarf+2)l"IPd(n;x;a) ^^^^^ 

{ad+iad+2, ad+iad+2/al; q)\n\ Hj^i {a'j+i/af, 
The main result in this subsection is the following theorem. 

Theorem 5.3. If{n,z,a) and {h,z,a) are related by (15. ip and if n,n G then 

Pd{n;x;a) ^ Pd{n;x;a). (5.5) 
Proof. Using (|2.ip and applying formula (j5.2p with fc, o, 6, c, d, e, / replaced by Uj, 

we can write the polynomial p„j in (|2.9p as follows 



V ^ ^9""^ZJ■+lg^"^^aJ+^,ao^"J9^^■-'^7^aJ+la7^^7^^J+l.„ ^ 

X 403 i-AT, -1 -1 TV, 1 -1 -1 ) 9' 9 



(5.6) 

Using the above formula for and (|5.ip one can check that the term of p„j 
in the variables (n, z,a) coincides with the 4^3 term of pfi^^j_^. in the variables 
(fi; z; ci) for j = 1, 2, . . . , d. Thus if we compute the ratio Pd{n; x; a)/Pd{h; x; a) all 
403 terms cancel and we get 

Pdin;i;a) ~ f}^ (l,)-".(ci,g^.-i5,,5,+ig^.-i5,+i/a2,a,+ia-i5,z-^_\;g)ft^ ' 

(5.7) 

Using (|5.ip we can eliminate Zj and Zj in the right-hand side of (|5.7p . From 
()5.1dP - (|5.1ep it follows that ad+iCk(i+2 = ced+ictd+2- Next, notice that ajZj = 
ad+iad+2q^'^'^'^~^ ■ Thus we have 

nfj _ TT \'^]'^d+l'^d+2) TT g "+1 3 J 
~-fij ^ 11 \nj 11 o-A'd+l-jMi ■ 

3 = 1 ^3 3 = 1 l"j"d+l°^d+2j J = l ^ 

It is easy to see that the second product on the right-hand side of the above formula 
is equal to one, by replacing Uj = Nj — Nj^i and fij = Nj — Nj^i. Therefore 

d -Mj d / -1 -1 \"j 



Using the identity 



3 = 1 ^3 ' 3 = 1 i"j"d+l"d+2j 



(«?';<?). = ^^^ for fc,ZGNo, 

(a;?); 



we obtain 



-pr (ajg^^-^Zj;g)n, _ -pr (ad+iad+2g^^-^+^''+^--' ; g)«, 

/J (<5,9^-^%;g)n, (ad+iad+2g^-^+^''+-^;g).., 

(ad+iad+2;g)jv,+iv,+i_, /J ("<i+i«d+2;9)jv,_i+jv,+i_. 
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It is easy to see now that the first product above is equal to 1, while the second 

(Qd+iad+2;ij)f 



simplifies to ("'^+i"''+^'^)"d _ Thus we have 



Similar manipulations show that 

d 



-pr {ajQ'^'-'zfq)n, ^ {ad+iad+2;q)\n\ 
^\{ajq^^-^zf,q)ii^ iad+iad+2; q)\A\' 



n 



_ {ad+iad+2/oil;q)\n\ -[4 ("l+i/aj ; g)n 



(5.10) 



n 



{ad+iad+2/al;q)\n\ ^J^ (<5j+i/a|; g)n,- ' 

The proof follows from equations ([SJ)) . (|5?8)) . (|5J)l and ((5?T0)) . □ 

5.2. The commutative algebra An- We denote by I?" the associative subal- 
gebra of I?^ of difference operators with coefficients depending rationally on the 
parameters aj . Clearly, the commutative algebra Az defined in Proposition 14.51 
is contained in P". Similarly, we denote by I?" the associative algebra of dif- 
ference operators in the variables n = (ni, 7^2, ... , Ud) with coefficients depending 
rationally on g"^ , g"^, . . . , g""* and the parameters aj. We use En^ to denote the 
forward shift in the variable n/c, i.e. for every function /(n) = . . . we 

have Enkf{n) = f{n + efe). 

Replacing j by d + 1 — fc in (|5.1ap we see that 

Otk+lZk+1 

From this equation it follows that a forward g-shift in the variable Zk will correspond 
to a forward shift in hd+i-k and a backward shift in hd+2-k for fc = 2, 3, . . . , d and 
to a forward shift in fid when k = 1. Thus, in view of the duality established in 
Theorem 15.31 we define a function b as follows 

b{ao) = ao (5.11a) 

b{a,) = ^oad+iad+2 ^,/2 for J = 1, 2, . . . , d + 1 (5.11b) 

O^d+2-j 

b{ad+2) = — g-1/2 (5.11c) 

b(z,) ^ "'^+2-3 ^»i+»2+-+«d+i-,-i/2 for J = 1,2,..., d (5.11d) 
aa 

b{E,.z^)^En,^,_^E~l^_^, forj = l,2,...,d (5.11e) 
with the convention that En^^-^ is the identity operator. 

Lemma 5.4. The mapping (j5.1ip extends to an isomorphism b : I?" — > I?", /n 
particular, the operators defined by 

2l^b{2l) for k = 1,2,..., d (5.12) 
commute with each other and therefore 

X = b(A)-C[£l\£^\...,£2] (5.13) 
is a commutative subalgebra ofD^- 
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Proof. Using ()5.1ip one can check that 

b(i?,,,J-b(/(z))-b(/(z<z^'=))b(i?,,.J, 

holds for every /c = 1, 2, . . . , d, which shows that b is a well defined homomorphism 
from T>" to V". The fact that b is one-to-one and onto follows easily. □ 

Recall that the eigenvalues Hj of the operators i2| are given in (j4.11bp . We 
denote 



(5.14) 



N /, ad+iad+2 \ /, oid+iad+2Zd+i-j 
Kj = b i^ij) = - 1 1 - 

\ Q.d+1-jZd+l-j J \ Oid+l-j 

T "d+i"d+2 , '2.ad+iad+2 

— — i 2 1 Xd+l-j- 

"d+l-J OLd+l-j 

With these notations we can formulate the main result. 

Theorem 5.5. The polynomials Pd{n;x;a) defined by equations (j2.9p and ([57 
diagonalize the algebras Az and An- More precisely, the following spectral equations 
hold 

£,^Pd{n; x; a) = ^jPdin; x; a) (5.15a) 
£,^Pd{n;x;a) = KjPd{n;x;a), (5.15b) 
for j — 1,2, ... ,d where £|, fij are given in (I4.1ip and Kj are given in (j5.12p . 

Remark 5.6 (Boundary conditions). Since £)! contains backward shift operators, 
and since Pd{n; x; a) is defined only for n € Nq it is natural to ask what happens 
when £^ produces a term with a negative Uj for some j. From (|5.11ep it follows 
that b{E^ 2) with i' E {—1,0, l}"* will contain a negative power of En^^^-j U — 2) 
in one of the following two cases: 

Case 1: Vj = 1, Vj-i = 0. In this case, h{E'^ ^) will contain E^^^^^,. Notice that 
the coefficient of El! , is Cy which has the factor ( 1 "^^^ — ) (in - see 

formula (|4?2b| ). Since = (1 - q-""^^-^-^) we see that this term is 

when nd_|_2-j = 0. 

Case 2: Vj — 1, Vj-i = —1. This time b(£'J' ^) contains E~^^^_,. The coefficient 
a has the factor = (l - ^-^71777) (l - S^T^) and KM ). 

Since b{Bjy^) = (1 — g^"<i+2-j)(l — (jfi-"d+2-3) ^fJ■Q gee that this coefficient is when 
nd+2-j = or 1. 

Proof of Theorem 15.51 Equation (j5.15a[) follows immediately from Proposition 14.51 
and the fact that Pd{n; x; a) and Pd{n; x; a) differ by a factor independent of z. It 
remains to prove (|5.15bp . Let us fix n,a. By Lemma for z £ (C*)'^ we can find 
(n, z, a) such that equations (|5.ip hold. If n € Ng then we can use Theorem 15.31 
and replace Pd{n;x;a) by Pd{n;x;a). Equation (I5.15bp follows from the fact that 
the operator £" in the variables n with parameters aj coincides with the operator 
£| in the variables z with parameters aj. We can use analytic continuation and 
complete the proof as follows. Fix n, a and write z in terms of the dual variables 
n, i.e. we put 

^ "rf+l«rf+2 gn,+n. + ...+n,^,_, fo^, fc = 1, 2, . . . , d. 
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Both sides of equation (I5.15bp become Laurent polynomials in the variables q^^ , 
(j"^, . . . , (j""*. Since (j5.15b|) is true for every ri € Ng, we conclude that it must be 
true for arbitrary n E C^, or equivalently, for arbitrary z e (C*)''. □ 
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